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Abstract. The theoretical interpretation of measurements of "wavefunctions" 
and spectra in electromagnetic cavities excited by antennas is considered. 
Assuming that the characteristic wavelength of the field inside the cavity is 
much larger than the radius of the antenna, we describe antennas as "point- 
like perturbations". This approach strongly simplifies the problem reducing the 
whole information on the antenna to four effective constants. In the framework 
of this approach we overcame the divergency of series of the phenomenological 
scattering theory and justify assumptions lying at the heart of "wavefunction 
measurements" . This selfconsistent approach allowed us to go beyond the one-pole 
approximation, in particular, to treat the experiments with degenerated states. 
The central idea of the approach is to introduce "renormalized" Green function, 
which contains the information on boundary reflections and has no singularity 
inside the cavity. 
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1. Introduction 

Cylindrical microwave cavities with different cross-sections and small height are now 
widely used as the simplest realization of "quantum" (wave) billiards [2 El [3] . Indeed 
for the lowest TM-mode (transversal magnetic) inside of the cavity the transversal 
component of electric field satisfies the two-dimensional Helmholz equation with 
Dirichlet boundary condition [4]. This problem is mathematically equivalent to the 
problem for the quantum particle inside a box with hard walls. This boundary problem 
can not be solved analytically for arbitrary cavity shapes, where the corresponding 
classical problem is non-integrable and shows chaotic behavior. Microwave cavities 
are a simple realization of "quantum" chaotic systems. In spite of its experimental 
simplicity, microwave billiards exhibit general properties of quantum chaotic systems 

Experimentally the cavity is always open, since one has to excite and measure 
the field inside. Usually one antenna is used for an excitation and the same or a 
second one for the measurement [5] (see FigfT]). These antennas can be considered as 
"probes" , assuming that the influence of antennas on the investigated properties is 
comparatively small. However this assumption quite often is violated. To understand 
the data in these cases we have developed a self-consisted theory of the experiment. 

The detailed analysis of the antenna construction represents an individual 
complicated problem [BJ. Therefore one has to treat antennas phenomenologically 
to obtain analytical formulas. The standard scattering matrix formalism [7J [HI [S] 
applied to the billiard system does not take properly into account singularities from 
the sources and contains divergent series. This disadvantage can be overcome if one 
uses an approach based on "point perturbations" of the system [9j [10] . The main idea 
of the approach is the introduction of the "renormalized" Green function instead of 
the Green function with coincident space arguments. 

In Section [2] we explain the idea of a point-like perturbation in two-dimensional 
case and derive the matching condition between wavefunctions inside of the antenna 
and inside of the cavity. Doing this we basically follow the line of Exner and Seba [11] . 
Ibidem we discuss the physical meaning of the coupling coefficients. In Section [3] we 
generalize the scattering matrix (^-matrix) obtained in [llj to the case of multiple 
antennas and estimate coupling coefficients for metallic and dielectric scatterers. 
Section[4]represents a general theory of an experimental wavefunction extraction based 
on the expression for ^-matrix obtained in the previous section. The theory is based 
on the one-pole approximation in the isolated resonances regime. At the end of the 
Section we show what happens when the "wavefunction" is measured in the case of 
degenerated states. 

While for closed cavities one can take into account only a finite number of 
resonances inside of the experimental range of wavenumbers, for open cavities this 
approximation fails. In Section [5] we show that for an open cavity with broad 
resonances the transmission coefficient between two antennas is proportional to the 
Green function of the unperturbed system. This gives the possibility to investigate the 
spatial structure of the Green function of an open cavity, for example nodal domain 
distributions jT2"llT3]. 

In Section [6] we consider the spectrum of a cavity with antennas connected by 
loops. Such a cavity represents a simple model of a quantum graph containing parts of 
different dimensions. In Section[7Jwe generalize the model to treat multimode cavities. 

|Appendix A| gives the link between the presented approach and the usual WW + 




Figure 1. a) The scheme of the physical cavity with antenna, b) Mathematical 
model of the cavity. 



approach where point-like perturbations are treated as (^-functions. In | Appendix B| 
we write the Green function for the cavity with antennas. Poles of this Green function 
determine resonances of the system disturbed by the antennas. 

2. Matching for point-like antenna 

We now discuss the general idea of the approach used in Ref. [11] which allows us 
to treat the presence of scatterers and antennas in a simple way assuming that the 
wavelength of the field is much larger than the characteristic size of the area of the 
attachment. Figure [1] illustrates the principal setup to measure spectra of a cavity. 
One sees that a quasi-one dimensional waveguide is attached to a cavity at some 
point. We will assume that the waveguide is straight and one can separate variables 
inside of it and define a scalar longitudinal wavefunction depending on the variable 
inside the waveguide. Furthermore we assume that the height of the cavity is small 
and the distance between transversal wavenumbers is large. Thus it is natural to 
suppose that the one-dimensional longitudinal wavefunction in the waveguide relates 
to the longitudinal wavefunction in the cavity in such a way that only some general 
and energy independent properties of the matching are important. Since the detailed 
construction of the matching usually is quite complicated, the only restriction applied 
is that the final operator describing the coupling between antenna and cavity should 
be linear and Hermitian. 

We substantially follow paper [IT] . Basically the idea of the approach is to 
construct a self-adjoint extension of the Hamiltonian in the presence of point-like 
perturbations. In various situations self-adjoint extensions were considered [101 114] 
and related problems were studied [El [16] . 

Let us consider the matching of a quasi one-dimensional antenna to a cavity. 
We denote by (pi(x), — oo < x < the longitudinal wavefunction inside the one- 
dimensional antenna and by ^(r), r = {x\,X2} a function inside the cavity. Thus 
ip = {ipi(x), (p2(r)} is a function on the space consisting of two parts: cavity and 
antenna. We assume that the antenna is coupled at the end x — to the point r = R 
of the cavity. We define the operator H which acts on the function <p in the following 
way: 

where <p"(x) = d 2 ipi(x)/dx 2 and A is the two-dimensional Laplacian. Since for TEM- 
modc in the waveguide and the lowest TM-mode in the cavity thresholds do not exist 
we did not consider them in Eq. ([T]). The generalization to higher transversal modes 
in the cavity will be presented in the Section [7] 
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In the stationary case we look for the solution of the spectral problem 

Hip = k 2 ip, (2) 

where k is the total wavenumber. 

Now we turn to the question how to define operator H to be Hermitian. This 
is equivalent to the introduction of proper boundary matching conditions for the 
functions <fii(x) and ^(r). 

In order to avoid divergent integrals at the origin first we consider the cavity with 
a small pricked ring of radius e with a center at the point r = R. We call this area 
f2 c . Following the general procedure [TTj, |H1 ED] we consider two quadratic forms 

(<p + Hi/>) =- f dx<pxi)'{ - [ d 2 r<p 2 A^ 2 , 

{^+H V ) t = - [ dxtw'i- I d 2 ri, 2 A V2 , (3) 

where bar means complex conjugation, superscript '+' denotes hermitian conjugation 
and angle brackets mean integration. Then consider the difference 

f-0 



J — OO 

+ / V{il> 2 V{p 2 - >p 2 Vtp 2 )d 2 r = [ipiipi - (piip'i\ x =o 



n,- 

dd[ih{d<p 2 /dp) - <f 2 (dtp 2 /d P )] p=e . (4) 



In the last equality we introduced cylindrical coordinates p,fi : {x\ = pcos$,x 2 = 
ps'm-d}, p = |r — R|, and assumed 

[ipi<Pi - <piip[] x =-oo = 0, 

[i> 2 v n cp 2 - ^ 2 v„V2] CT = o, (5) 

where a is the physical boundary of the cavity and V„ is the derivative along the outer 
normal to the boundary. This is equivalent to the assumption that the functions or 
their normal derivatives are zero at the boundary (Dirichlet or Neumann boundary 
conditions or their linear combination). 

Assuming that functions <p 2 , ip 2 have the same kind of singularities for r — > R as 
the Green function for the two-dimensional plane [TH [19] 

G(r,R;fe 2 ) = ^(fc|r-R|), 

G(r - R R; fc 2 ) - -i- (in + 7 ) + j, (6) 

where is the Hankel function of the first kind and 7 ~ 0.5772 is Euler constant, 
we can write asymptotic expansions for p — ► 0: 

Va ( P , <p) = (2tt)- 1 l^ Hp) + Lf + f v (p, V ), 

Mp, <p) = {^r x Lt Hp) + Lf 2 + Ufa <p), (7) 

where Lq 2 , Lf 2 , Lq 2 , Lf 2 are constants. We have written the factor (27r) _1 in order 
to simplify the further calculations. Here p is calculated in arbitrary length units. 
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Though ln(p) depends on the choice of the length unit, the sum of the first two terms 
in expansions ([7]) is invariant. We will discuss this feature later. 

Assuming that f v ^{p) — > and \df v ^,/dp\ < const for p — > 0, we get 

r2ir 

L/ H)lp=e — A 

Finally in the limit e — > equality (|4]) takes the form 



lime / d0[M9<h/9p) ~ ^(d^/dp)] p=e = L$*L? - L* 2 L^.(8) 



(ip + Hilj) - (ip+H<p) = lira - (^ + Hip) f 

= - + Lf 2 Lt 2 ~ LfL?. (9) 

Let us introduce vectors r] v = (^(0), Lf 2 ), = (<Pi(0), L% 2 ), r]^ = (ip[ (0), Lp), 
Cv> = (^i(O)j Lq 2 ) and rewrite the last equality in the form 



(V + Hi>) - {iP+H v ) = {n^Ct) - {Cp,rn,) , (10) 

where (•, •) denotes the scalar product of two vectors. Since the operator H is linear 
we write 

r\ v = H v , m = A C-0 (11) 

and get that matrix A has to be Hcrmitian. Thus the general matching condition 
reads 



(12) 



^(0) \ _ (A B\( ^(0) 
Lf )-\C D ){ LI 2 

where A, D are real and C = B. As it is assumed that H is fc-independent, the 
constants A, B, C and D are fc-independent as well. From the last formula one can 
see that for B = C = antenna and cavity are decoupled, which means there is no 
mixing between functions ipi, (f2- fi(0) = A(fi(0), Lf 2 — DLq 2 . Thus we conclude 
that B, C are coupling constants, A describes the boundary condition at the end of the 
antenna and D describes the scattering strength of the antenna. A — oo corresponds 
to <pi(0) = (Dirichlet boundary condition) and A — corresponds to (f'i(0) = 
(Neumann boundary condition). 

Let us now prove the invariance on the length scale of the expression (J7| . From 
the matching condition (fT2|) we have 

L\ 2 = DI% 2 +CVi(0). (13) 
Substituting this relation in (JT]) we get 

<p 2 (p, <p) = Lt 2 [(27T)- 1 ln(p) +D}+ CVi(0) + f v (p, <p). (14) 
Let us introduce the scattering length (3 by the following equality 

D = -^-lnp. (15) 

Then (fT4")) takes the form 

T<P2 

<p*0>, <p) = in j 3 + dpi (o) + Up, <p), (16) 

which is invariant. The definition of the scattering length (|15|) corresponds to the 
condition that the wavefunction is zero at the distance /3 from the perturbation point 
when B = C = 0. It differs from the definition given in [5]. 
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3. 5-matrix for a cavity with multiple antennas 

Here we consider the case with v antennas attached to a cavity at points 
Ri, R.2, . . . , R„. We look for the solution in the form 

<pu(x) = a l e lkx + b t e- zkx , , , 

^(r)=Er=iWr,R j; fc 2 ), [U) 

where ifu (x) is the solution inside the i-th antenna, ip2 (r) is the field inside the cavity, 
G(r,Ri; k) is the Green function in the cavity and Fi are unknown coefficients which 
have to be defined. 

The Green function G(r,Ri; k 2 ) has singularities at the points r = R^. Thus in 
the neighborhood of R^ we write 

G(r - R,, R, l5 k 2 ) - --L In |r - R,| + £(R i5 k 2 ), (18) 

Z7T 

where £(Ri;fc 2 ) represents the renormalized Green function. It plays the key role 
in the following calculations. Numerically the renormalized Green function can be 
calculated by means of the boundary element method [20 . 

From the definition we see that £(Ri; k 2 ) depends on the chosen system of units. 
However all equations determining physical observable quantities will of course be 
invariant. 

Then we define generalized boundary values Lq? , L\l at the matching points 



<^ 2i (r -» Ri) -> (2n)- L L^ ln(|r - R 4 |) + L 



li ' 



rV>2 _ 77 



Lif? =F i e(R l ;fc)+J]^G(R l ,R,;fc 2 ). (19) 

Using matching condition (fT2|) . we get 

ik(a,i -bi) \ _ ( Ai Bi \ ( ai + b 

F^(R i ;fc 2 ) + E^^G(R 4 ,R J ;fc 2 ) )~\Oi D t J\ -Ft 



(20) 



These equations can be rewritten in a matrix form 

BF = (A - ik)a + (A + ik)b, (21) 

{G + D)F = C(a + b), (22) 

where vectors F, a, b have elements Fi, ai, bi respectively, matrix G has elements 

Gij =G y (R 1 ,...,R,;fc 2 ) = | ^ k 2 ), l^'j, ( 23 ) 

and diagonal matrices A, B, C, D have elements Ai, Bi, d, Di on diagonals 
respectively. Substituting F from (|2"2")l into (|21[) we get 

(M - ik)a + (M + ik)b = 0, (24) 
where M = A— B(G + D^C. Thus 

b = Sa, (25) 
where the S'-matrix has the form 

M + ik y ' 
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The Green function of a closed cavity is real, thus from Eq. 
S-matrix is unitary. Using the equality 

M -ik A-ik _ . 1 1 



2ik- 



M + ik A + ik 
we represent S in the following way: 

A-ik 2ikB 



M + ik 



(M - A) 



A 



one sees that the 



(27) 



S = -- 



A + ik A + ik 



G + D- 



BC 

A + ik 



C 



A + ik 



(28) 



For the case of a single antenna formula (|2"5|) was obtained in [TT] . 

From (|28|) we see that poles of S'-matrix k = k(R\, . . . , R„) are the complex 
solutions of the equation 

G(Ri,...,R,;fc 2 ) + J D h< 



det 



0. 



(29) 



A + ik 

In the limit of a closed system (B, C — > 0) we get the equation determining the 
spectrum of the system with only scatterers inside 

r 2N 



det[G(Ri,...,R„;Ar) + D] = 0. (30) 

Eq. ([3171) shows that in the limit D — > oo, poles of the system approach the poles of 
the unperturbed system. 



5. i. One- dimensional cavities 

If an antenna is attached to a narrow waveguide with a width comparable to the 
wavelength, then the cavity can be considered as the one-dimensional one. The S'- 
matrix for such a cavity can be calculated analogously to two-dimensional case. We 
introduce the operator inside the cavity 

-cp'Kx) 

(-cP/dxl+kl)^^) 



Hip = 



(31) 



where fcf is a threshold from the quantization in x 2 direction. We are looking for 
the solution of the spectral problem @. We denote X the position of the antenna 
attachment, Xi, X r are the coordinates of the ends of the cavity, X\ < X r . Calculating 
the difference of two quadratic forms we get 



(ip + Hili) - (ip+Hip) 



dx{-(p x ^'{ + '4) l (p1) + 



x- 



X,. 



dxi(-ip2ip2 + VWa) + / dxx{~ip 2 ip2 + faf'i) = 



X+c 



where e — * 0, 



^ = (pi(0),^- £ ) 1 ^ + £ )) 1 

= (M0),Mx - e)^' 2 (X + e)), 
^ = (^(0),^(X-e)> 2 (X + e)). 



x- 



X+e 



(33) 
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Then we write i)u>^ = AC,^. From (|32[) A is Hermitian. Usually one assumes that 
tp(X - e) = tp(X + e) = <p(X), thus 




A = | C D I ) . C-S4) 
This gives 

^ 1 (0)=A^ 1 (0) + S^ 2 (X), (35) 

^ (X - e) - tft 2 (X + e) = Ctpx (0) + D^ 2 (X). (36) 

For a cavity with a scatterer inside (B = C = 0) we get 

^ (X - e) - <p' 2 (X + e) = Jty a (X) . (37) 

Exactly the same condition holds for the derivatives of the solution of the following 
equation with true delta-scatterers 

[-d 2 /dx 2 - D6(x - X)]<p 2 = (k 2 - kl)<p 2 . (38) 

We have mapped the problem of a ID cavity with a point-like scatterer to a Schrodinger 
equation with a ^-potential. Such a mapping is possible in one-dimensional system 
only. 

From the asymptotic 

V 2 (x 2 ^X)^L% 2 \x 2 -X\+L^, (39) 

the condition (|37|) and the assumption that Lq 2 (3 + Lf 2 = we find D — 2/(3, where 
(3 is a scattering length. 

Let us now construct the S'-matrix for a number of antennas attached to a cavity 
at the points X±, X 2 , ■ ■ ■ ,X V , Analogously to (|T7|) . (TSU|) we write 

<pu(x) = a l e lkx + b t e- zkx , . . 

V 2(x 2 )^j::=iF l G(x 2 ,X l ;k 2 ), (W) 

where k\ = k 2 — k 2 and 

<p' u (0) = Ai<pn(p) + Bi<p 2i (Xi), 

V2 (X i -e)-<p 2 (X i + e)=C i <p li {0) + D i <p 2 (Xi), [ ' 

where G(x 2 , X; k\) is the Green function of the closed one dimensional cavity. Using 
the identity 

x-t 

= 1, (42) 



lim 



-^G(x 2 ,X;k 2 1 ) 
dx 2 



X+e 



we find analogously to 



A-ik 2ikB 



A + ik 



C 



A + ik' 



(43) 



A + ik A + ik 
where 

G ij =G{X i ,X j ;4). (44) 

On the contrary to ([3H| poles of the closed system (B, C — > 0) tend to the 
unperturbed ones when D — > 0. 
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3.2. Example: Single antenna attached to an infinite cavity 

One can imagine a Gedankexperiment where an antenna is attached to an infinitely 
large cavity. In such a cavity there are no reflections from the walls. Experimentally 
this situation can be realized if one places a perfect absorber along the boundary 
of the cavity. In terms of the cross-section such a resonator can be considered as a 
two-dimensional plane. 

Let us consider a single antenna attached to a plane at the point R. In this case 
the unique element of S'-matrix 5*11 is equal to the reflection coefficient and A, B, C, D 
are scalar constants. From ([28} we find 

A-ik 2ikBC 
11 ~ ~ A + ik + (A + ik) 2 

where 

£(R;fc 2 ) = -i-(ln|+ 7 )+l. (46) 

The last equality follows from ©. In the limit of a closed cavity the pole of Su is 

k = 2ie~~*/f3~ 1.123i//3, (47) 
where (3 is the scattering length. The eigenvalue of the problem |(3J) is k 2 = 
—Ae~ 2l (3~ 2 ~ — 1.261/3~ 2 . This eigenvalue corresponds to the eigenstate induced by a 
point-like scatterer with a wavenumber below the threshold. This means that a point- 
like scatterer in two-dimensional case can be always treated "quantum-mechanically" 
as an attractive well. Since in electromagnetic cavities the threshold for the lowest 
transversal mode is equal to zero, this pole can be observed only for higher transversal 
modes. 



f(R; k 2 ) + D 



BC 



A + ik 



(45) 



3. 3. Estimation of the scattering length 

In the previous example we showed that a point-like perturbation can be always treated 
as an attractive well. Nevertheless using it one can model not only dielectric but also 
metallic scatterers. The apparent contradiction is explained by the fact that in the 
long-wave limit the whole information about the scatterer is reduced to its scattering 
length which can be the same for the attractive as well as for the repulsive scatterer. 
This description fails whenever two scatterers are close to each other or close to the 
boundary (compare [H]). Below we estimate a scattering length for attractive and 
repulsive scatterers. 

Let us consider a cylindrical cavity with radius R and a cylindrically symmetric 
scatterer of radius ro whose center coincides with the centrum of the cavity. We shall 
treat the scatterer in a quantum-mechanical way as a potential and shall consider 
states with zero angular momentum. 

We put the origin of the frame in the center of the cavity. We denote 
the wavefunctions inside and outside of the scatterer by ip in (r;k) and i/j out (r; k), 
respectively. Here r = |r|. Then 

Vw(r; k) = ci(fc) Jo(fcr) + c 2 (k)N (kr), (48) 
where Jo and Nq are Bessel and Neumann functions of zero order. At the boundary 
of the cavity ijj ou t{R;k) = 0. This gives c{k) = c 2 (fc)/ci(fc) = -J Q (kR)/N (kR). 
Equating logarithmic derivatives at the boundary of the scatterer we find 

ntVi^T-Qifc) _ -fcr [Ji(fcr ) +c(k)N 1 (kr Q )} 
tpin(ro;k) J (kr ) + c(k)N (kr ) 
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where ip'. = dipm/dr. We assume that kr <C 1. Then one can put k — on the left 
hand-side of P5| and use the asymptotic expansions Jq{x) — 1— x 2 /A+0(x 4 ), N (x) = 
(2/7r)[ln(a;/2)+7] + 0(a; 2 ln(a;)), J^x) =x/2 + 0(x 3 ), N^x) = -2/(ttx) +0(xln(x)) 
jS) . This gives 

r o ^ ra (r o ;0) c(fc) 

lMm;0) 7r/2 + c(fc)Kfcr /2)+ 7 ]' 1 j 

From the last equation we get the condition 

C(fc) = - ln(fcr /2) - 7 + Oli 0)/M4>o; 0)] ' (51) 
Let us consider what happens if we simulate the scatterer by a point-like 
perturbation. Using the following equalities 

G(r, 0; k 2 ) = - 1 [J (fcr)/c(fc) + 2V (fcr)] , 

e(0;fc 2 ) = -T47T-^^|+7V (52) 



4c(fc) 2tt V 2 

we find from Eqs. (HJ), (J30j) 



This gives the following quantization condition 



«*> = -«T? <54 » 



Comparing Eqs. (|51[) and (|5l)l we find 



^ = ^oxp(-^#4,). (55) 



nr^- n ( r o;°), 

Let us discuss the last result. First it shows that a singular point-like scatterer can 
model an attractive as well as a metallic scatterer. For the metallic scatterer of the 
radius ro we find ipi n (ro) — and f3 — ro in accordance with (|16p . 

Let us consider now an "attractive scatterer" . For simplicity we consider a 
potential well of constant depth — uq. Then ipi n (r;k) = aJo(nr), >c = (uq + fc 2 ) 1 / 2 , 
ipin(ro;0)/[roip' in (ro;Q)] = -Jo(v^n))/[v^ r o^i(V*to r o)]- If ^fmr Q < 1 then 

fe(^o;Q) 2_ 

ro^„(r o ;0)- M0 r 2 ' 1 j 

Substituting (f56|) in ([55]) we find that the scattering length f3 can increase exponentially 
for an attractive scatterer. For quite small ro we have kf3 3> 1. Thus we conclude 
that exponential increase of the scattering length for an attractive scatterer leads to 
the new characteristic length different from the radius of the scatterer. 

From (f5"3"]) we find that the cases of small "metallic" scatterers (k(3 <§C 1) and 
small attractive ("dielectric") scatterers (k(3 ^> 1) lead to different signs of the shift 
of resonances of the resonator with respect to the empty one. For metallic scatterers 
resonances arc shifted to higher wavenumbers and for the dielectric ones they are 
shifted to lower wavenumbers (see also [5]), which will become evident from the next 
section. 
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Figure 2. Experimentally measured resonances for the rectangular cavity 
(47.3 X 20.0 cm) with metallic and teflon scatterers. The vertical dotted line 
corresponds to the eigenwavenumber of the rectangle. The solid curve shows 
the modulus of the reflection coefficient from the empty cavity with a single 
exciting antenna. The dashed and the dashed-dotted curves show the resonances 
in the presence of a cylindrical metallic (rrj = 1.50 mm) and dielectric (teflon, 
ro = 2.56 mm) scatterers positioned at x\ = 33.4 mm, X2 = 50 mm, respectively. 



4. Wavefunction measurements 



In this section we consider the technique of wavefunctions measurements in the closed 
cavity. If the cavity is closed, its spectrum is discrete. The presence of antennas opens 
the cavity, thus all eigenvalues become resonances of the iS-matrix. If some resonances 
are isolated from the rest, they can be considered separately and the influence of all 
others can be neglected. The same consideration takes place for an open cavity with a 
number of sharp resonances, isolated from the rest. The wavefunction corresponding 
to an isolated sharp resonance can be measured experimentally. Below we consider 
the idea of wavefunction measurements and some effects appearing if resonances are 
overlapping. 



4-1- One pole approximation 

Close to an eigenvalue of the discrete spectrum the Green function may be written 
as 

GM;^),^», (57) 
JSp - k z 

provided r and R are far enough from each other. Function £(R; k 2 ) can be written 
in the form 

e(R;fc 2 ) = -|^ + e(R;fc 2 ), (58) 

where £(R; k 2 ) has no pole at k 2 = E^, In one pole approximation (|58|) reads 

ibf. (R) 

^R;k 2 )^^-^+a^E,). (59) 
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Using (T5T)) , ([55)1 one can calculate the inverse matrix appearing in Eq. (|2"5)) [2] . Indeed 
in this case 

G + D-^^ = a\^)(^\+V, (60) 

where a = (E^ — k 2 )~ 1 , = ^(R,), V = D - BC/{A + ik), D VJ = A% with the 
renormalized dimensionless coefficient 

D i = D i + £(TL i ;E fl ). (61) 

Using the identity 

a\^)(ij\+V l + a^V- 1 ^)' K ' 

we can rewrite Eq. (f2"5]l for k 2 close to in the form 



s = _A-ik li | 2ikf,B 



c 



(63) 



S M -A; 2 + {^|y- 1 |V). 
where fc M = \JE~^. This immediately gives us poles of the iS-matrix: 

*-B r -±- ^ . (64) 

Vl A-l-B,p/(A, + *„) 

Since nondiagonal elements SV,-, i 7^ j are proportional to ipiR^ipCRj), they contain 
an additional information on the sign of the wavefunction. 

Formula (|64|) implies the possibility to measure wave junctions through the shifts 
of resonances at different points of the cavity induced by a point like perturbation. In 
Ref. pQ the author used a movable point-like scatterer and fixed antenna whereas in 
[H [22] the antenna has been moved. 

Figure [2] illustrates shifts of the resonance of the empty rectangular cavity (solid 
curve) as well as in the presence of a small metallic (dashed curve) and a dielectric 
(dashed-dotted curve). One sees that for a metallic scatterer the resonance is shifted 
to a higher wavenumbeiQ while for the dielectric scatterer it is shifted to a lower 
wavenumber. This observation is in agreement with the arguments of the Section [3. 31 



4-2. m-poles approximation 

In m-poles approximation one takes into account only contributions from m poles in 
some range of wavenumbers. For r / R we get an approximate expression for the 
Green function truncating the convergent series at the mth term: 

GM;^EW (65) 

£(R; k 2 ) can be written in the form 

£(R;fc 2 ) = £-p^ + £(R;fc 2 ), (66) 

% For the larger metallic scatterer (ro = 2.28 mm) we observed a shift of the resonances to the lower 
wavenumbers in the low-fc regime (k < 0.9cm" 1 ). This might be related to the fact that in the 
lower wavenumbers regime the distance between antenna and scatterer becomes comparable to the 
wavelength and the scattering length of the antenna. Then the proper description of the situation 
requires to go beyond one-pole approximation or even beyond the concept of a point-like scatterer. 
This experiment has to be investigated further in details. 
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where £(R; k 2 ) has no poles for E\ > k 2 < E m . The approximation assumes that the 
energy dependence of £(R; k 2 ) can be neglected. Then we get approximately 

C(R; k 2 ) ~ + «R; S), (67) 

jU=l M 

where E is some value from the range E\< E < E m . Using ([551) . (E3 we can wr he 
G + D = W + {H - k 2 Y x W + D. (68) 

Here H^, = £^<W> W^i = W (**■»')' D %j = A = A + £(R;;£)- We see that 

inside of the given window of wavenumbers ^(Rjj-B) is an offset renormalizing the 
scattering strength. 

Let us write S-matrix (|28|) in m-poles approximation. We neglect energy 
dependence everywhere apart from the Green function. We have 



G + D- B ° 



1 



A + ik E 

where k E = y/E, V = D- BC/{A + ik E ). We write 

1 

V + W + (H - k 2 )~ 1 W ~ 



W+(H -k 2 )-^ + V 



(69) 



(70) 



Finally we get 



V~ 1 -V~ 1 W + WV^ 1 (71) 

H-k 2 + WV- 1 W+ { ' 



A — ik E BCV- 

S = ; h llk E - 



A + ik E *(A + ik E ) 2 

2ik E BV- 1 , 1 V~ X C , . 

W + W 72) 

A + ik E H - k 2 + WV~ 1 W+ A + ikY 

_i _ D( A 2 + k 2 E ) -BC(A + ik E ) _ 

(AD - BC) 2 + k 2 E D 2 



y-l = ^ = 1 _ X ( } 



If D > 1 then 



^-D- 1 , VT 1 ~ - kB ° . (74) 

R 1 D 2 (A 2 + k 2 E ) V ; 

Introducing an effective Hamiltonian 

H cS = H + WV^ 1 W + -iWV I - 1 W + , (75) 

we can write 

A- ik B BCV- 1 2ik E BV- 1 TI , 1 Tzr V^C 
S =—a T^ + ^Tl ^7 t— — r, W+- -W- — . 76 

A + ik E (A + ik E ) 2 A + ik E H eS -k 2 A + ik E 

Experimentally one usually removes the global phase of S'-matrix, describing the 
antenna matching itself. Theoretically this means that instead of (|76p one introduces 
a new matrix S' 

A + tk Ec< 1 o ., BCV- 1 2ik E BV-\ ir+ 1 jjr V~ X C 

s =--a — = l-2ik E -^ — pT+ * ~, — w t? t^ w ~a — i— ■ i 77 ) 

A~ik E A 2 +k% A- ik E H eS - k 2 A + ik E 

Poles of S'-matrix can be found from the equation 

det[fZeff- k 2 ] = 0. (78) 
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For the closed system the last equation reduces to 
det[£T + WL>- 1 W + - k 2 ] = 0. 



(79) 



Effective Hamiltonians of the type (|75|) were used in [3 [3 [8] . Resonances of S- matrix 
in a general situation are discussed in [23] . 

4-3. Experimental "wavefunctions" in the case of degenerated states 

Let us consider a cavity with the degenerated eigenvalue = E 1 /i = 1, . . . , to. In 
this case any weak perturbation mixes all eigenstates corresponding to E. In the 
wavefunction measurement an antenna is a perturbation. Moving the antenna one 
changes mixing coefficients and produces a new field configuration inside the cavity. 
Therefore it is not clear whether one can easily treat the experimental wavefunction 
in this case. We show below that the experimental pattern in such a situation is just a 
sum of squares of all eigenfunctions ip^ (r) . This combination is invariant with respect 
to the choice of eigenfunctions. 

The simplest approximation in this situation is to take into account only m 
eigenfunctions corresponding to E. Let us calculate an experimental pattern which 
can be reconstructed from the measurement with a single antenna. The resonance 
position for the fixed antenna position R can be calculated from Eq. ([75)1 . Using 
Eq. (J67J we get 



From Eq. ([5Tj) we conclude that in the case of degenerated eigenvalue the 
experimentally measured pattern should correspond to the sum of squares of 
corresponding eigenfunctions. 

In Fig. [3j experimental examples from a rectangular microwave billiard with a 
side ratio of 2:1 are presented. For the technical details the reader is referred to 
Ref. [24j . The figure shows experimental and theoretical patterns corresponding to 
doublets consisting of the 15th, 16th states and the 25th, 26th states. States of each 
doublet are degenerated, the data were obtained from the reflection measurement with 
a single movable antenna. Thus one expects to measure the sum of squares of two 
eigenfunctions. Pictures a)-d) illustrate experimental data, all others show theoretical 
simulations. 

Fig. [3^) shows the resonances corresponding to the 15th state with the antenna 
at the positions marked by a square, a triangle and a disk at Fig. [3J;). Solid, dashed 
and dashed-dotted lines in Fig. [3^) correspond to a square, a triangle and a disk, 
respectively. Here only a single resonance is observed and no splitting could be 
seen. Fig. \3]p) shows the shift of the resonance while the antenna is shifted along the 
horizontal line in Fig. [5J;). In Figs. [3J:) and d) the obtained resonance wavenumber 
shifts patterns are plotted in a gray-scale. In Figs. (3^)-h) the squares IV'il 2 , l^l 2 
of the rectangular eigenstates are plotted. In Figs. G3)-j) \^pi\ 2 + |V>2| 2 is plotted, 
demonstrating the perfect agreement with the experimental results shown in Figs.[3j;)- 




(80) 




(81) 



d). 
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Figure 3. Comparison of experiment and theory for two degenerated states of 
the rectangular billiard with side ratio 2:1: a) experimental resonances for the 
degenerated 15th and 16th states for three antenna positions marked by a square, 
a triangle and a disk in c); b) the resonance wavenumber corresponding to 15th, 
16th states versus antenna position on the horizontal line in c); c) resonance 
wavenumber-shift pattern in a gray-scale; e) l^ll 2 ; g) I'fel 2 ; i) l^il 2 + IV^I 2 for 
the 15th and 16th states; d), f), h), j) the same for the degenerated 25th and 26th 
states 



5. Transmission measurements in open systems 

Recent measurements of nodal domains in open microwave cavities [121 113] raised 
a question on the theoretical description of this type of experiments. The crucial 
difference compared to the case considered above of a closed cavity consists in the 
absence of sharp resonances of the S'-matrix. This makes it impossible to use the one- 
pole approximation being the base of wavefunctions measurements in closed systems. 

In the experiment the transmission coefficient from the fixed antenna to the 
movable one is measured and interpreted directly as a "wavefunction" corresponding 
to the continuous spectrum of the unperturbed open system. We show below that in 
fact in the experiment the Green function of the cavity is measured. To prove this 
statement we start from Eq. (|28[). The straightforward calculation gives for the case 
of two antennas 

^(R.R^) = - x(Al+ f k ^ A2+ik f (^k% (82) 
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A = A(Ri,R 2 ;fc 2 ) = dct 



G(Ri,R 2 ;fc 2 ) + L> 



BC 

A + ik 



Let us assume that k is far from all poles fc M and scattering is weak, i.e. 

Di = Di- (2^)- 1 [ln(fc/2) +7] = -(27r)- 1 [ln(A:/3 J /2) + 7] 
Here we used (|15p . Thus for quite distant antennas we can write 



(83) 



(84) 



(85) 



Measuring Si 2 (Ri, R2; k 2 ) at a fixed wavenumber and a fixed antenna position R2 
we get a value proportional to G(Ri, R2; k 2 ). An additional phase shift in S\% due to 
the prefactor of the Green function in Eq. (|82|) can be easily removed for the closed 
cavity, since the Green function in this case is real. 



5.1. Example: Single antenna attached to a semi-infinite cavity 

The following example illustrates general effects taking place in the vicinity of a 
boundary. We consider two antennas attached to a semi-infinite cavity at points 
Ri, R 2 . The cross-section of the cavity is a half plane. 

Let the boundary of the cavity coincides with a;2-axis in the plane. We define 
Ri = (#1, x 2 ), R2 = (ii, x 2 ). The transmission coefficient S\2 through these antennas 
is defined by (|82|) where G(Ri, R2; k 2 ), £(Ri;fc 2 ) can be found from the method of 
images: 

i 



1-H, 



(!)/ 



G(R 1 ,R 2 ;fc 2 ) = -^(fc|R 1 -R 2 |)-i/^(fc|R 1 -R 2 |) , R 2 = (-x u x 2 

£(R i; fc 2 ) 

For simplicity we assume fc|R! 
In this approximation 



(86) 



1 

'2tt 



k \ i 



(2fcri)J. (87) 
R 2 | 3> 1 and neglect nondiagonal elements in G in 



B1C1 



C(R 2 ;fc 2 ) + ^ 2 - 



B2C2 



Ai+ikJ y y ' * A 2 +ik / 

If the movable antenna is placed far from the boundary then kx\ 3> 1. Using the 
asymptotic expansion of the Hankel function 10] and (|15|) we write 

e(Ri; k 2 ) +D X ~-— ln(A/3i) + - (l 7== )■ (89) 



2ir 4 V ^/nkxi ' 

Thus we see that at distances of the order of the wavelength from the boundary A 
contains an oscillating term. This term decays when one moves the antenna away 
from the boundary. The origin of this oscillating part is the mirror source of the 
movable antenna. This mirror source can be treated as an "impurity" charge, thus 
the oscillating behavior of A may be interpreted in terms of Friedel oscillations [53] . 

Let us now approach the movable antenna to the boundary. If kx\ <S^. 1 then 
using ©, (fT5|) we write 

^(R 1 ;fc 2 ) + ^ 1 ~^fln^+ 7 V (90) 



2vr 



The comparison of the last formula with an assumption (|85p shows that at distances 
less than the wavelength the measured quantity can differ from the Green function. 





Figure 5. Shift of the same resonance of the rectangular cavity as shown in Fig. [2] 
due to a loop connecting two antenna at positions (14.0, 8.0) cm, (40.4, 5.0) cm. 
Solid line: both antennas are closed by 50 f2 load at the end, dashed line: the 
antennas are connected by a wire. Dotted line: eigenwavenumber of the rectangle. 



6. Cavity with closed loops 

In this section we consider a cavity with closed loops (see Figure [4}. This model is 
another realization of a quantum graph 26J. In order to simplify calculations we do 
not consider attached antennas. 



6.1. Cavity with short loops 

Analogously to Eq. ^ we introduce the Hamiltonian inside the cavity 

Hip = -Aip. (91) 
Repeating steps (|5|)- (fTU)) wc get 

(tp,HiJ>) - WW) = (%,Cv) - (C v ,%) , (92) 
where r) v , ( v , r^, Cv> are vectors with components rf^ = = Lfg, rj 1 ^ = 

rV>2 /■ _ T 4>2 
^li ) S<A — ■ 

Analogously to we write 

V v = D ( v , Vl> = D U- (93) 
Here D is a Hermitian matrix with elements Dy, i, j = l,...,u, where v is the 
number of inserts including both scatterers and points of loops attachments. Elements 
Dij = 0, i j if points i and j are not connected. In contrast to (f!?5|) matrix D 
contains non-diagonal elements describing loops. 
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Using matching condition (|93p we can find eigenwavenumbers of the closed cavity. 
We look for the solution in the form 

V 

^(r)=^^G(r,R i; fc 2 ). (94) 

i=l 

Substituting ([19]) in ([93]) we get 

[G(R x ,...,TL„;k 2 )+D]F = 0, (95) 
where F — (Fx, . . . , F v ). Thus the equation for the eigenvalues fc 2 is 

det[G(Ri,...,R„;fc 2 )+D] =0. (96) 

6. 2. A more realistic model of the loop 

The model of loops considered above is appropriate only for short loops or for some 
narrow energy region since it does not contain the spectrum of the one-dimensional 
wire connecting different points. To take into account the "inner structure" of the 
loop let us consider a more realistic model. We consider a cavity with a loop of the 
length L matched to it at two points Ri, R2. We neglect the curvature of the wire. 
Analogously to Eq. (fT7]) we write 

ifii(x) = ae lkx + be~ lkx , 

^(r)=E?=i FiGir^k 2 ). 

Matching conditions (|12[) at points Ri, R2 give 

ik(a -b) \ _ f A x B x \ / a x + b x 

F x ^B. x ;k 2 )+F 2 G(K x ,K 2 ;k 2 ) ) ~ \ C x D x ) \ -F x 

ik{ae tkL - be- lkL ) \ _ ( A 2 B 2 \ ( ae lkL + be~ lkL 

F 2 ^(R 2 ;fc 2 )+F 1 G(R 2 ,R 1 ;fc 2 ) ) ~ \ C 2 D 2 ) \ -F 2 

From Eq. ([9"5]) we find that the spectrum of the system is defined by the equation 

det (I 1 2 , „ ) = 0, (99) 



(97) 



(98) 



v Z 2 G + D 
where Z x , Z 2 are 2x2 matrices 



, Ai-ik A x +ik , 

{A 2 -iky kL (A 2 +ik)e-* kL ( " l01 



(101) 



C x C x 
C 2 e lkL C 2 e-' LkL 

If the cavity and the wire are decoupled then B = 0, Z 2 = and Eq. ([M]l reduces to 

det(Zi)det(G + L>) = 0. (102) 

The last equation gives the unification of spectra of the cavity with scatterers inside 
and the wire. 

In Fig. O one sees the experimentally observed shift of the same resonance of 
the rectangular cavity as shown in Fig. [2] The dotted line shows the position of 
the eigenwavenumber of the rectangular cavity. The solid line corresponds to the 
modulus of the reflection coefficient from the cavity excited by the antenna with two 
additionally attached antennas each closed by a 50 load. The dashed line illustrates 
the modulus of the reflection coefficient from the same cavity, were the antennas were 
connected by a wire. 
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7. The S'-matrix for multimode cavity 



If there are more than one open transversal modes in the cavity with transversal 
wavenumbers < k, n — 1, ...,iV the ansatz of the solution (jTTJ) has to be 

generalized. To treat this situation let us consider N identical "copies" of the cavity 
(see figure[6]) with wavefunctions ^(r), ■ ■ ■ ( PA r +i( r ) and introduce the operator acting 
on the spacial variables of the antenna and the N cavities: 



(-A + (fc«) 2 V 2 (r) 
V (-A+(fef>)>iv +1 (r) / 



(103) 



where k" is n-th transversal wavenumber. We look for the solution of the spectral 
problem ([2]) in the form 

ip u (x) = a t e lkx + he~ lkx , 

^ +1 (r)=^^G(r,R i; (fc^) 2 ), 



where k 



(n) 



[k 2 



(fci n) ) 2 ] 1/2 , n 



1. 



(104) 

, N (for the lowest TM-mode in 



electromagnetic billiards fc^ 1 ' = and kiP = k). 

The same procedure as for a single transversal channel gives us the following 
expression for the ^-matrix: 

M - ik 



S = 



M + ik' 
where M = A — Z, Aij = AiSij, 



(«) m-ii( nn ')^»( n ') 



(105) 
(106) 



Here the matrix Q has elements 

n (nn') _ A(n) c 
Win ~ ^ii u nn' 



D {nnl) 8- 



(107) 



where G^ 



6y(R,,R i; (*j n) ) a ) 
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8. Conclusions 



In the paper we presented the general selfconsistent theory for measurements on 
electromagnetic cavities with point-like antennas and scatterers. In the framework 
of the considered approach we justified the idea lying on the heart of wavefunctions 
measurement. The experimental patterns for the cases of degenerated states are 
explained. The generalization of the basic idea allowed to describe cavities with more 
than one open transversal channel and antennas connected by loops. Mathematical 
objects appearing in the last cases can be considered as quantum graphs. 
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Appendix A. Operator algebra for point scatterers in a closed cavity. 

Similar to the one-dimensional case one usually treats a point-like scatterer as a 5- 
potential which can be written in an arbitrary representation as a matrix WW + , 
where W is a coupling matrix [5]. However in a multidimensional case a point-like 
perturbation does not correspond exactly to a ^-potential. We show below that instead 
of a 5-potential in the two-dimensional case a point-like scatterer situated at the point 
R is described by a product S(r — ~R)1Z, where 1Z is a renormalizing operator. 
From Eqs. (J7J), (fT9)) we find that the function ip 2 obeys the equation 

{H-k 2 )ip2=^2Fi6(r-Ki). (A.l) 
»=i 

Using Eq. (12]) with B = C = we write 

(H - k 2 )ip 2 {r) = -J2 D^5(r - R;)7^ 2 (r) (A.2) 

i=l 

where 

nip 2 {v) = lim [(1 - ln(|r' - r|)(r' - r)V)(/3 2 (r')]- (A.3) 

r' — >r 

Thus 1Z acts on the function <p(r) without a singularity at the point r as an identical 
operator and gives the renormalized value for a function with a logarithmic singularity 
(see Eq. (0). 

In the basis of eigenfunctions Vv( r ) ^ ne matrix elements of the S- function are 
given by 

(lMr)|*(r - R)|Vy(r)) = ^(R)^(R). (A.4) 
In this basis Eq. (|A.2[) reads 

(H- fc 2 )|^ 2 > = -WD- l W + Tl\^), (A.5) 
where \ f2)^ = (Vvl 1 ^), W^i = ^(Rj). The solvability condition of Eq. (|A.2[) is 



det 



H - k 2 + WD- 1 W + K 



= 0. (A.6) 
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The last equality can be rewritten in the form 



det 



H 



WD-i-W+ll 



det [if - k 2 } x 



x det 



1 + WD- x W + Tl{H - fc 2 )- 1 



det [if - k 2 } x 



x det[D _1 ] det 



D + W + K{H - k 2 )-^ 



0. 



(A.7) 



For any vectors \<p), we have (ip\lZ\<p) — (ip\<p) if (ip\<p) exists. Thus for i ^ j 
we have 



(A.8) 

(A.9) 
(A.10) 



[w+n(H k 2 r x w] t] = y MgiMiM , 

From this we conclude that 

W + K(H -k 2 )-^ = G. 
Finally the equation for the spectrum reads 

det[fi-/c 2 ]det[G + L>] =0. 
The last equation is in accordance with (f5TT|) . 

Appendix B. Green function for the system with antennas. 

Let us consider a cavity with a number of fixed antennas. In order to measure 
the wavefunction of such a system one has to introduce a movable antenna. The 
wavefunction of the cavity with a movable antenna is the Green function of the cavity 
with only fixed antennas. Then the natural question is to calculate Green function of 
the cavity with a number of fixed antennas. 

Let us look for the Green function of the system with fixed antennas G p in the 

form 



G p (r,R; k 2 ) = G(r, R; k 2 ) +^F J G(r,R :( ;fc 2 ), 



(B.l) 



where G(r, R; fc 2 ) is the Green function for the cavity without antennas. The function 
(|B.1I) satisfies the equation for the Green function. The Green function should 
correspond to outgoing waves spreading inside of the fixed antennas, thus inside of 
the ith fixed antenna the solution is 

<p u (x) = b t e- lkx . (B.2) 

In the neighborhood of R^ the Green function has the following asymptotic 

G p (r -> Ri , R; fc 2 ) -> G(R l; R; fc 2 ) + 



+ Y / F 3 G(R t ,Ry,k 2 ) + F l 

Using the condition ([12")) we get 

— ikbi 

G(R,,R; fc 2 ) + E j¥i R,.R,: k 2 ) 



-^lndr-R.D+^fc 2 ; 



At Bi 



-Fi 



(B.3) 



(B.4) 
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From the last equality analogously to (|21[) . (f2"2"| we find 



G p (r,R; k 2 ) 



G(r,R; k 2 ) - ^ G(r, R,; k 2 ) 



G + L> 



5G 



A + ifc 



G(Rj,R; fc 2 ).(B.5) 



If B = C = then the last formula reduces to the so-called "Krein's resolvent formula" 
for the Green function of the system with point-like scatterers [10] . Krein's formula can 
be also obtained by the operator technic of Appendix 1 using the following identities 

G p (r,R;fc 2 ) ' 



(Rl 



(Rl 



1 



1 



-W 



H - k 2 + WD~ 1 W+TZ I ' 

1 „.^„ 1 



H-k 2 H-k 2 W + TZ(H — k 2 )~ 1 W + D 



W + K- 



H-k 2 



(B.6) 



where |r) M = tp^(r), |R) M = ?/v( R )- 
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